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O ; ABSTRACT 

. We present new constraints on the evolution of dark energy from an analysis of Cos- 

mic Microwave Background, supernova and X-ray galaxy cluster data. Our analysis 
I/"") , employs a minimum of priors and exploits the complementary nature of these data 

sets. We examine a series of dark energy models with up to three free parameters: the 
current dark energy equation of state wq, the early time equation of state w e t and the 
scale factor at transition, at. From a combined analysis of all three data sets, assuming 



> 

in 



Of 



a constant equation of state and that the Universe is flat, we measure wq = — 1.05io 12- 
Including w ct as a free parameter and allowing the transition scale factor to vary over 
. the range 0.5 < at < 0.95 where the data sets have discriminating power, we measure 

f*^) ' wq = — 1.27ig'3g an d w ct — — 0-66ig'g2- We find no significant evidence for evolution in 

, the dark energy equation of state parameter with redshift. Marginal hints of evolution 

C^) ■ in the supernovae data become less significant when the cluster constraints are also 

included in the analysis. The complementary nature of the data sets leads to a tight 
constraint on the mean matter density, Q m and alleviates a number of other parameter 
degeneracies, including that between the scalar spectral index n s , the physical baryon 
density fl^h 2 and the optical depth r. This complementary nature also allows us to ex- 
amine models in which we drop the prior on the curvature. For non-flat models with a 
constant equation of state, we measure wq — —1.09~t [-£ and obtain a tight constraint 
on the current dark energy density, fide = 0.70 ± 0.03. For dark energy models other 
than a cosmological constant, energy-momentum conservation requires the inclusion 
of spatial perturbations in the dark energy component. Our analysis includes such 
perturbations, assuming a sound speed = 1 in the dark energy fluid as expected for 
■ Quintessence scenarios. For our most general dark energy model, not including such 

perturbations would lead to spurious constraints on w c t which would be tighter than 
those mentioned above by approximately a factor two with the current data. 

Key words: cosmology:observations - cosmology:theory - cosmic microwave back- 
ground - supernovae - x-ray clusters - dark energy 



1 INTRODUCTION 

The precise measurement of the Cosmic Microwave 
Background (CMB) made with the Wilkinson Mi- 
crowave Anisotrop y Probe (WMAP) llHinshaw et al-llioO^: 
iKoeut et al.ll200al has improved our knowledge of a wide 
range of cosmological parameters. However, a number of de- 
generacies between parameters exist which cannot be broken 
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with current CMB data alone and which require the intro- 
duction of other, complementary data sets. 

Some of the most important parameters and degenera- 
cies concern dark energy and its equation of state. Since 
observations of distant type la supernovae (SNIa) first in- 
dicated that the expansion of the Univer se is accelerating 
iRiess et alJl!99d iPerlmutter et al]|l999lh there has been 
enormous interest in this topic. The most straightforward 
way to incorporate an accelerated expansion into cosmolog- 
ical models is by adding a constant term to the Einstein 
equations - the cosmological constant. However, this leads 
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to an extreme fine tuning problem wherein one must ad- 
just the initial density of this constant to 10 _120 Mp 1 in nat- 
ural Planck units. To alleviate this, a scalar field model, 
dubbed Quintessence, was introduced which, when the po- 
tential is ca refully chosen, can avoid the fine tuning of initial 
conditions JPeebles fc Ratral Il98& iRatra fc Peebles! Jl98St 



Wettericbll988l:lFerreira fc Jovcell 998: Caldwell ct al. 1998; 



Zlatev et al.lll999T) . When describing the background evolu- 
tion of the Universe with such models, it is sufficient to 
know the equation of state for the dark energy i.e. the ratio 
of pressure and energy density, w = pde/pde- Whilst a cos- 
mological constant has w = — 1 at all times, for most dark 
energy models the equation of state parameter is an evolving 
function of redshift, w = w(z). 

In order to learn more about the origin of cosmic accel- 
eration and dark energy, it is crucial to constrain the evolu- 
tion of the dark energy equation of state. In the first case, 
this requires us to examine whether the accelerated expan- 
sion can be described by a cosmological constant or if there 
is need to go beyond this description. Most early attempts 
to parameterise the evolution of dark energy were carried 
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redshift SNIa and the arrival of other, complementary con- 
straints means th at we can now start to ask the same ques- 
tions of real data <Knop et aljl2003l : IRiess et al.ll2004) . 

The data for SNIa can be used to measure the luminos- 
ity distances to these sources independent of their redshifts. 
This constrains a combination of the dark matter and dark 
energy densities in a different way to observations of CMB 
anisotropies. The combination of the two data sets is there- 
fore useful in breaking parameter degeneracies. However, 
the simplest, linear expansion in redsh ift for the dark en - 
ergy equation of s tate advocated by e . g.lMaor et ail ||200 it) . 
lAstierl feOOll) and I Weller fc Albrechd d 200ll I 2002^ cannot 
be applied to the high redshifts probed by the CMB. For 
this reason, the linear para mete risation was extended by 
IChevallier fc Polarskil (1200 ll) and lLinderl <l2003fi to a model 
in which the equation of state at low redshifts (late times) 
wo, and at high redshifts (early times) w ct , could be speci- 
fied separately. Although more suitable than the low redshift 
linear expansion for the analysis of CMB data, this param- 
eterisation also has a short-coming in that the transition 
between wo and w e t always occurs at redshift z — 1 and 
with a fixed transition rate, which is not representative of 
the full rang e of scalar field dark energ y models of inter- 
est (s ee e. gJWeUe^^^Albrechd ll 20021) '). ICorasaniti et all 
(2003) and lBassett et all i2004f) extended this prescription 
further, allow ing the transition to o ccur at an arbitrary 
time and rate. lCorasaniti et alJ ll2004t) applied this extended 
model to a combination of SNIa, CMB and galaxy redshift 
survey data, noting the presence of strong degeneracies be- 
t ween a number o f the derived parameters. The analysis 
of ICorasaniti et ail ll2004f) included a limited exploration of 
models with an equatio n of state w < —1: so called phantom 
models llCaldwel]ll2002ri . While the extension to w < — 1 is 
challe nging in terms of the physics involved (ICarroll et alJ 
2003), it is interesting from a phenomenological point of 
view, particularly given that the best-fit to current super- 



novae data is obtained for models with w < — 1 iRiess et alJ 
12004 . 

It is important to note that the analysis of 
ICorasaniti et al] (12004]) did not include dark energy pertur- 
bations for models crossing w = — 1. While a cosmological 
constant is spatially homogeneous, this is not true for an ar- 
bitrary dark energy fluid or Quintessence. One must include 
perturbations in the dark energy component, not just for 
consistency reasons but also because the exclusion of them 
can lead to erroneously tight constraints on w from large- 
scale CMB anisotropies dWeller fc Lewisll2q03h . 

It was realised by Maor^t alJ f2001f) and 



I Weller fc Albrechd (I 200 ill 2002T) that in order to constrain 
the evolution of the equation of state with supernovae 
observations, it is necessary to use a tight prior on the mean 
matter density of the Universe, f2 m . Recent measurements 
of the gas fraction in X-ray luminous, dynamically relaxed 
clusters made with the Chandra X-ray Observatory provid e 
one of our best constraints on fl m dAllen et al] 12004]) . 
These data also provide a direct and independent method 
by which to measure the acceleration of the Universe, 
providing additional discriminating power for dark energy 
studies. As we shall demonstrate here, the combination of 
CMB and X-ray cluster data can also play an important 
role in breakin g othe r key parameter degeneracies (see 
also lAllen et al.1 ]2003ll 'l. For these reasons, we have used a 
combination of X-ray gas fraction, CMB and SNIa data in 
this study. 

In the following sections we first introduce our choice of 
parameterisations for the dark energy equation of state. We 
then discuss the individual data sets and how they probe 
cosmology. Our results are presented in Section [I] Sectional 
discusses the results and summarises our conclusions. 



2 DARK ENERGY MODEL 

A number of different parameterisations for the evolution 
of the dark energy equation of state parameter, w(z), have 
been discussed in the literature. The si mplest is the lin- 
ear parameterisation: w(z) = wp + w ' z iMaor et al]l200ll : 
IWeller fc Albrechdl 200lL I 20021 lAstierl 1200 it) . However, as 
mentioned above, this model is not compatible with CMB 
data since it div erges at high re dshift. IChevallier fc Polarskd 
d200l]) (see also iLinderl J2003)) proposed an extended pa- 
rameterisation which avoids this problem, with w(z) = 
Wo + w\z/{l + z). This model can in principle be used 
to distinguish a cosmological constant from other forms of 
dark energy with a varying w. However, this parameterisa- 
tion is not repre s entati ve of standard Quintessence models. 
ICorasaniti et al] dl2003l) proposed a generalised parameteri- 
sation which is better suited to the problem. However, that 
model includes four parameters and exhibits large degenera- 
cies when applied current data. 

Here, we use an exten sion of the model discussed by 
IChevallier fc Polarskil |200J) and lLinderl <2003ft ■ which stops 
short of the full extension suggested bv ICorasaniti et al] 
(2003). The primary short-coming of the parame t erisat ion 
proposed bv IChevallier fc Polarskil <200ll) : iLhiderj ll2003ft is 
that it uses a fixed redshift, z = 1, for the transition between 
the current value of the equation of state and the value at 
early times, w ct = Wo + Our model introduces one extra 
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parameter, z t , the transition redshift between w c t and wq, 
such that 

w e tZ + wozt w ct (1 — a) at + Wo(l — a t )a 

W= : = i ,/ - ; '— , (1) 



z + z t 



- a t 



a(l — 2a t 

where at is the transition sca le factor. (The parameterisa- 
tion of ICorasaniti et, all l)200ot) also introduces an arbitrary 
transition rate between wo and w s t-) 

Energy conservation of the dark energy fluid results in 
evolution of the energy density with the scale factor, such 
that 



Pdc(a) = Pdo,oa 3 e 3 ^ 



(2) 



where pd c ,o is the ener gy density of the dark energy fluid 
today. Using the parameterisation of equation Q we obtain 



f W ^ a / ^ da — w ct In a + (w ct — wo)g{a; at) , (3) 
Ji a 

a(l - a t ) 



with 



a(l-2a t ) 



- at 



(4) 



Setti ng Zt = 1 or a t = 1/2, we recover the parameteri- 
sation of iLindeil i2003T) . Hereafter, we shall refer to this as 
the zt ~ 1 dark energy model. From the Friedmann equa- 
tion, the evolution of the Hubble parameter H(z) — HqE(z) 
is given by 



E(z) = Vfim(l + + fidc/(*) + «k(l + Z) 2 , (5) 

with 

f(z) = (1 + 2 )3(l+™ct) e -3( I ,„t-»o)9( Z ;^) ) (g) 

where fi m , fide, fik are the matter, dark energy and curva- 
ture densities in units of the critical density. 

3 DATA ANALYSIS 

We have performed a likelihood analysis using three cos- 
mological data sets: CMB, SNla and the X-ray cluster gas 
fraction. 

For the CMB ana lysis we have modified the camb 1 
code llLewis et a l. 2000) to include the relevant dark energy 
equation of state parameters. For the calculation of CMB 
spectra, we have accounted for the effects perturb ations in 
the dark energy component JWeller fc Lewisll2003l) . We as- 
sume that the sound speed of the dark energy fluid, c 2 = 1, 
a choice t hat is well motivated for standard Quintessence 
scenarios dWeller fc Lewisll2003l) . although some other well- 
motivated dark energy models such as k-essence scenarios 
iArmendariz-Picon et all l2000t) include an evolving sound 
speed. We note the presence of an extra term in the per- 
turbation equations due to the variation of the equation of 
state with time, which sources the density perturbation with 
the velocity perturbation . This effect will be dis cussed in a 
forthcoming publication jRapetti fc WelleJl200j). 

We use three CMB data sets: WMAP llverde et alJ 
l2003l:lHinshaw et all2003tlKogut et"ai]l2003h (including the 
temperature-polarization cross -correlation d ata) , the Cos- 
mic Background Imager (CBI) dPearson et al.l2003T) and the 

1 http://camb.info 
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Figure 1. The 68.3 and 95.4 per cent confidence limits in the 
(n m ,tuo) plane for the various pairs of data sets and for all three 
data sets combined. A constant dark energy equation of state 
parameter is assumed. 



Arcminute Cosm ology Bolometer Array Receiver (ACBAR) 
JKuq et al.ll2004h . The latter data sets provide important 
information on smaller scales (I > 800). 

For th e SNl a analysis, we use the gold sample of 
iRiess et al.l (j2004fl . marginalising analytically over the ab- 
solute magnitude M as a "nuisance parameter". We fit 
the extinction-corrected distance moduli, fio — m — M = 
51ogd_L + 25, where m is the apparent magnitude and cLl is 
the luminosity distance in units of Mpc defined as 



d,L = tt- "Lfj sinh 



Hoy/ilk 



where fi^ 



fin 



fide 



dz 



(7) 



For the X-ray gas mass fraction analysis, we use the 
data and method of lAllen et alJ (120041) . fitting the apparent 



redshift evolution of the cluster gas fraction with the model 



6 fi b 



l+O.wVh) fi 



di CDM (z) 



df(z 



(8) 



where djf(z) and d|| CDM (z) are the angular diameter dis- 
tances (dA = dz,/(l + z) 2 ) to the clusters for a given dark 
energy (de) model and the reference standard cold dark 
matter cosmology, respectively, fib is the mean baryonic 
matter density of the Universe in units of the critical den- 
sity, Ho — 100 ftkmsec -1 Mpc -1 and & is a bias factor 
that accounts for the (relatively small amount of) baryonic 
materi al expelled f r om g alaxy clusters as they form. Fol- 
lowing lAllen et alJ i2004) . we adopt a Gaussian prior on 
b — 0.824 ± 0.089, which is appropriate for clusters of the 
masses studied here. Note that the prior on b includes a 10 
per cent allowance for systematic uncertainties in the nor- 
malisation of the /gas (2) curve, although we note that even 
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Figure 2. (Left panel) The marginalised constraints on wo assuming a constant dark energy equation of state. The solid line is for all 
the three data sets together, the short dashed line for clusters+CMB, the long-dashed line for SNIa+CMB, and the dot-dashed line for 
clustcrs+SNIa (with priors from HST and BBN). The right panel shows the marginalised constraints on Q m for each combination of 
data sets. 



doubling this s ystematic uncerta inty has only a small effect 
on the results lAllen et alJl2004) . 

We have included our extension of the CAMB code into 
the COSMOMC package, which provides an efficient sampling 
of the posterior likelihoods using a Markov Cha in Monte 
Carlo (MCMC) a lgorithm 2 |Lewis fe Bridlel2002l) . We have 
als o included th e]RJ^s_ et al, " ^(XMlT siiDernovae sample and 
the lAllen et alJ (120041) gas fraction data into the analysis, us- 
ing the COSMOMC code to calculate the posterior probability 
densities. 

For our standard analysis we have varied nine "cos- 
mological" parameters: the baryon density Qbh 2 , the cold 
dark matter density Qdmh 2 , the Hubble constant Ho, the 
reionisation redshift z rc , the spectral index n s , the ampli- 
tude of the fluctuations A s and the dark energy equation 
of state parameters Wo, w e t and at. We assume zero ten- 
sor components and a negligible neutrino mass. The bias 
parameter 6 associated with the X-ray cluster data is an 
additional parameter in the fits. We have marginalised ana- 
lytically over the intrinsic magnitude M of the supernova*. 
Except where stated otherwise, our analysis assumes that 
the Universe is flat (f2k = 0). For the analysis of the clus- 
ter data without the CMB data, we use Gaussian priors 
on n b h 2 = 0.0214 ± 0.0020 from Big Ba ng Nucleosynthesis 
(BBN) constraints llKirkman et all200&t and h = 0.72±0.08 
from o bservations made with the Hubble Space Telescope 
(HST) iFreedman et al.ll200lf) . 

For each different model+data combination, we have 
sampled four independent MCMC chains. The length of 
these chains vary from the simplest models of constant equa- 
tion of state with around six thousand samples per chain up 
to our most general model with around one hundred and fifty 
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Table 1. The median values and 68.3 per cent confidence intervals 
from the analysis of the various pairs of data sets and for all three 
data sets combined, assuming a constant dark energy equation of 
state. The last column states the x 2 P er degree of freedom for 
each combination of data sets. 



Data combination 



wo 



X 2 /dof 



Cl+CMB 


-1 9 o+0.17 


U.AO4_ 034 


1467.4/1378 


C1+SN(+HST+BBN) 


-1 04+°' 14 


n 9ss+ - 043 


210.2/178 


SN+CMB 


_1 04+° 13 
LMq: -0.16 


U.aiZ_Q Q44 


1625.9/1508 


Cl+SN+CMB 


1 - uo -0.12 


n 2QC .+0.031 


1652.5/1534 



thousand samples per chain. We have ap plied the Gelman- 
Rubi n criterion to test for convergence jGelman fc Rubml 
1992). Convergence is assumed to be acceptable if the ra- 
tio of the between-chain and mean-chain variances satisfies 
R — 1 < 0.1. We have also checked for convergence by ensur- 
ing that consistent final results are obtained from sampling 
numerous small subsets of the chains. 



4 DARK ENERGY CONSTRAINTS 

We have employed a series of different parameterisations for 
the dark energy equation of state: (i) w constant (ii) a model 
with wo and w e t free, but with the transition redshift fixed 
at Zt = 1 (at = 0.5) (iii) a similar model with the transi- 
tion fixed at z t = 0.11 (a t = 0.9) (iv) a similar model with 
the transition fixed at z t = 0.35 (at = 0.74), which approx- 
imately splits the cluster and supernovae data about their 



© 2004 RAS, MNRAS 000.HHTT1 



Constraining dark energy with galaxy clusters, supernovae and the CMB 5 




Figure 4. The 68.3 and 95.4 per cent confidence limits in the (Qm',wo ,w e t) plane for all three data sets combined using various fixed 
values for the transition redshift. The solid lines show the results on (f2 m ,ii>o)- The dashed lines show the results on (O m ,ui e t)- The left 
panel is for z t = 1 (at = 0.5), the centre panel z t = 0.11 (a t = 0.9) and the right panel z t = 0.35 (at = 0.74). The uncertainty in w e t is 
much larger than for wo in the left panel, which reflects the paucity of data at high redshifts. The z% = 0.35 transition splits the cluster 
and SNIa data into similarly sized low and high redshift subsamples. The horizontal dotted line denotes the cosmological constant model 

(«0 = Wet = -1). 
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Figure 3. The 68.3 and 95.4 per cent confidence limits in the 
(fi m ,uio) plane for the various pairs of data sets and for all three 
data sets combined. The z t = 1 dark energy model is assumed. 

median redshifts and (v) a model in which the transition 
redshift is a free parameter. 

Figure shows the constraints on wo and f2 m for the 
constant dark energy equation of state model. We see that 
the combination of the three data sets leads to tight con- 
straints on Wo and fl m , which are in good agreement with 
the cosmological constant scenario (wo = —1). This figure 
also demonstrates the complementary nature of the con- 
straints provided by the various pairs of data sets, in par- 
ticular SNIa+CMB and clusters+CMB. 

The power of the combined clusters+SNIa+CMB data 
set is also evident in the marginalised probability distribu- 



tions for wo (left panel of Figure^ and Q m (right panel). 
The marginalised 68.3 per cent confidence limits on wo and 
S7 m for the various data combinations, assuming a constant 
dark energy equation of state, are summarised in Tabled 

Figure |3] shows the results obtained using the zt = 1 
dark energy model. Comparison of Figs and shows how 
the joint confidence contours on (O m ,wo) open up when w e t 
is introduced as an additional free parameter. This is partic- 
ularly prominent for the clusters+CMB combination where 
the wo region is extended into the positive branch due to 
the degeneracy between wo and w e t, discussed below. The re- 
sults on wo and U! e t for the zt = 1 model are shown in the left 
panel of Figure [I] The figure shows how, for this model, the 
present data constrain wo more tightly than w e t- This simply 
reflects the paucity of cluster and SNIa data at redshifts be- 
yond Zt = 1. Note that the CMB data provide an upper limit 
of Wet J; at high redshifts; for w ct > 0, the dark energy 
component would become significant at early times, causing 
modifications to the predicted CMB anisotropy spectrum. 

It is important to recognise that the choice of transition 
redshift, z t = 1, described above is arbitrary. We have there- 
fore examined the constraints obtained for other values of Zt 
(a t ). The centre panel of Figure 2] shows the results using a 
late transition model with z t = 0.11 (a t = 0.9). We see that 
the cosmological constant (wo = w ct = — 1) again lies within 
the allowed 68.3 per cent confidence (la) regions. Unsurpris- 
ingly, the constraints on w ct in the late transition case are 
better than for the Zt = 1 model, reflecting the presence of 
more cluster and SNIa data beyond the transition redshift. 
Naturally, this at the expense of a weaker constraint on wo- 
lf we select a transition redshift close to the median 
redshift for the SNIa and cluster data sets, one might ex- 
pect to obtain comparable constraints on wo and iw e t. In 
principle, this approach could provide improved sensitivity 
when searching for evolution in the equation of state pa- 
rameter. (In detail, we expect the constraints on wo to be 
slightly better than those for w e t using the median redshift 
model, since the precision of the individual cluster and su- 



© 2004 RAS, MNRAS OOO.miTTl 



6 D. Rapetti, S. W. Allen and J. Weller 



Table 2. The median parameter values and values at the peaks of the marginalised probability distributions (and 68.3 per cent confidence 
intervals) for the various dark energy paramctcrisations, using all three data sets combined. Results are listed for both flat and non-flat 
priors. The last column states the \ 2 P er degree of freedom for each parameterisation. 



equation of state 


median wq 


median w et 


median f2 m 


Peak too 


Peak w Bt 


Peak Q m 


X 2 /dof 


constant (flat) 
constant (non flat) 
z t = 1 (flat) 
z t = 1 (non flat) 
split, zt = 0.35 (flat) 
arbitrary zt (flat) 


i - uo -0.12 
-1 09+°- 12 

-1 io+°- 23 

± - lu -0.19 
-l.U8_ 23 

-i ,o+ - 34 

!.OU_q 2g 

i 97 +0.33 
± z -0.39 


-0.87±° : g 
1 90+0.86 

± - zo -2.06 
-0 61 +°- 40 

-0 66+ 43 


fl 295+ 031 
u.z»o_ 027 

o 14 +0.040 

300+ 029 
u.ouu_ 02g 

fl Q9Q+0.046 
U.328_ 04Q 

n o nn +0.028 

299+ 029 
u.zyy_ Q 027 


-1.04 ±0.10 

-1 02+ - 08 
± - uz -0.20 

-1 16+ 022 

_1 14+0.31 

_1 49+0-41 
1, * a — 0.20 

1 90+O.34 
lzJ -0.46 


-0.05«;? 7 
-0 09+ - 12 

uuy -2.16 

-0 m+ 09 

u - iu -0.71 

-0 12+ - 11 
u - lz -0.76 


0.293 ±0.028 
31+ 04 

0.30 ±0.03 

n 09+0.04 

0.30 ±0.03 
0.298 ±0.028 


1652.5/1534 
1651.3/1533 
1650.6/1533 
1648.1/1532 
1649.0/1533 
1648.9/1532 
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Figure 5. The 68.3 and 95.4 per cent confidence limits 
in the (wo,w e t) plane obtained using the SNIa±CMB (blue, 
thin solid), Clusters+CMB (green, long-dashed) and Clus- 
ters±SNIa(±BBN±HST) (magenta, dotted) data sets. The bold 
lines show the results for all three data sets combined and the 
dark circle marks the cosmological constant model. The dashed 
line shows the no evolution models (wo = Wet). The transition 
redshift is fixed to zt = 0.35 in all cases. 



pernova measurements are lower at high redshifts.) The right 
panel of Figure^Jshows the results obtained fixing z t = 0.35 
(at = 0.74), a redshift close to the median redshift for both 
the cluster and SNIa data sets. In this case the uncertainties 
on wo and to e t are indeed similar and the combined size of 
the confidence regions is reduced. However, the cosmological 
constant remains an acceptable description of the data. The 
marginalised results on wo, w ct and fi m are summarised in 
Tabled 

Within the context of searching for evolution in the dark 
energy equation of state, it is interesting to note the con- 
straints that arise from the combinations of SNIa±CMB, 
clusters±CMB and clusters±SNIa(±BBN±HST) data sep- 
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Figure 6. The distribution of MCMC samples in the (wo,w e t) 
plane for the case of a transition scale factor allowed to vary over 
the range < a t < 1. The colours indicate the value of at. 



arately. Figure shows the results in the (iBo,w e t) plane for 
the three pairs of data sets and for all three data sets com- 
bined, using the z t = 0.35 dark energy model. We see that 
the combination of SNIa+CMB data provides marginal evi- 
dence for evolution in the equation of state, in that the 68.3 
per cent confidence contours exclude the no evolution line 
(wo = Wet)- However, the clusters±CMB and clusters±SNIa 
data are consistent with the cosmological constant at the 
68.3 per cent level, and the effect of combining all three data 
sets (bold contours in Figure^ is to remove the marginal 
evidence for evolution hinted at in the SNIa±CMB data 
alone. 

Figure |S] clearly shows the degeneracies between wo and 
w c t for the different data combinations, and demonstrates 
the importance of including the CMB data when attempting 
to obtain the best constraints on wo and w c t- Comparing 
the dotted contours [clusters±SNIa±(BBN±HST)] with the 
bold contours [clusters±SNIa±CMB] one sees that as well as 
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Figure 7. (Left panel) The 68.3 and 95.4 per cent confidence limits in the (flm;wo,w c t) plane for a general dark energy model with the 
transition scale factor allowed to vary over the range 0.5 < at < 0.95. The solid lines show the results on (n m ,u>o). The dashed lines 
show the results on (n m ,ui c t). The horizontal dotted line denotes the cosmological constant model (u>o = w et = — 1). (Right panel) The 
constraints in the (wo,w e t) plane for the above model (bold, solid lines). For comparison, the constraints obtained for the various fixed 
transition redshifts are also shown: Zt = 1.0 (green, dashed), zt = 0.35 (blue, thin solid), Zt = 0.11 (red, dotted). The dark circle marks 
the cosmological constant model. 



providing a tight upper limit on w e t as discussed above, the 
inclusion of the CMB data also leads to tighter constraints 
on Wo and the exclusion of large negative values for w e t- 

The most general dark energy model we have examined 
includes wo, w c t and the transition scale factor, at, as free 
parameters. In the first case, a t was allowed to take any 
value in the range < a t < 1. The distribution of MCMC 
samples from this analysis is shown in Figure |5] This figure 
re-emphasises the point that when the transition in the dark 
energy equation of state occurs at late times (low redshifts; 
light coloured sample points) w ct is confined to a relatively 
narrow band (—2 < w ct < 0) and the constraint on wo 
is poor (—2.5 < wq < —0.5). When the transition occurs 
earlier (at high redshifts; darker points), the constraint on 
wo is improved and the constraint on w ct is weakened. In this 
case, some sample points even populate the region beyond 
w ct > 0. [If the transition occurs at sufficiently early times, 
even the CMB data cannot provide a tight limit on w ct . 
It is likely, however, that an early equation of state with 
w > 1/3 violates BBN bounds, if the dark energy component 
is significant at t he time when the BBN species freeze out 
jBean et alJ^OOll) .] 

It is clear from Fig|S|that interesting constraints on both 
wo and w e t can only be obtained when the transition redshift 
is restricted to lie within the range spanned by the cluster 
and SNIa data. Otherwise large peaks in the marginalised 
probability distributions will occur that will simply reflect 
an inability to distinguish between models with transition 
redshifts beyond the range of the present data. For this rea- 
son, we have carried out a second analysis in which a t was 
allowed to vary only over the range 0.5 < a t < 0.95; a sen- 
sible compromise given the current cluster and SNIa data. 
Fig|7|shows the confidence contours in the (n m ;wo,w e t) (left 
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Figure 8. The marginalised results on wo (solid line) and w e t 
(dashed line) for the general dark energy model with the transi- 
tion scale factor allowed to vary over the range 0.5 < at < 0.95 
(c.f. Fig|7J. 

panel) and (tflo,w e t) planes (right panel, bold-solid lines). 
The marginalised results on too and w e t are shown in Fig|H| 
Again, the results obtained with our most general dark en- 
ergy model are consistent with a cosmological constant. 

As mentioned above, our analysis accounts for the ef- 
fects of spatial fluctuations in the dark energy component. 
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Figure 9. The 68.3 and 95.4 per cent confidence limits in the 
(wo, w e t) plane obtained from analyses which account for (blue, 
solid contours; as in Fig|7J or incorrectly neglect (red, dashed con- 
tours) the effects of dark energy perturbations. The model used 
incorporates a free transition redshift for the dark energy equation 
of state and is fitted to all three data sets: clusters+SNIa+CMB. 

Most previous studies have not accounted for these per- 
turbations, despite the fact that this leads to violation of 
energy-momentum conservati on whenever dark energy is 
not a cosmological constant llHul I2004I: ICaldwell fc Doranl 
I2005I) . It is important to ask whether the inclusion of these 
perturbations has a significant effect on the results; it has 
been argued by some authors that dark energy perturbations 
can be neglected if the equation of state remains around 
the cosmological constant value. However, we find that for 
an evolving equation of state, neglecting the effects of such 
perturbations can lead to spuriously tight constraints on the 
dark energy parameters. 

For a constant equation of state. IWeller fc Lewisl <|2003T) 
showed that the inclusion of dark energy perturbations 
leads to an opening up of the (fi m ,«;o) contours, allowing 
more negative values of u>o- Repeating their analysis us- 
ing our three data sets (clusters+SNIa+CMB) we measure 
a reduced effect, due to the complementary nature of our 
data sets. Neglecting the effects of dark energy perturba- 
tions leads to only a small shift in the marginalized prob- 
ability distribution for wq and slightly tighter constraints 
(wo = — 0.991q'^[; see Table Q for the results obtained in- 
cluding perturbations) . 

For our most general, evolving dark energy model, how- 
ever, the effects of perturbations in the dark energy compo- 
nent are more important and neglecting them can lead to 
spuriously tight constraints. Figure [5] compares the results 
in the (wo,w t) plane obtained when including dark energy 
perturbations (solid contours) or neglecting them (dashed 
contours). When the effects of dark energy perturbations 
are wrongly ignored, we obtain spuriously tight constraints 
on w = — 1. 25^^34 an d especially w ct = — Q-OTtj^J; the 
apparent uncertainties in the latter are reduced by a factor 
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Figure 10. The 68.3 and 95.4 per cent confidence limits in the 
(f2 m ,u>o) plane assuming a constant dark energy equation of state 
and co mbining CMB data either with the 2dF galaxy redshift 
survey iPercival et alj[2002l) (red, solid contours) or with the X- 
ray cluster gas mass fraction data (green, dashed contours). 

of ~ 2 from the values in Table [5] Figure [5] also shows that 
when we incorrectly neglect the effects of dark energy per- 
turbations, the sharp boundary provided by the CMB data 
set around Wet ~ is not reached. 

Finally it is interesting to compare the constraints on 
the matter content and dark energy parameters from this 
study with other work where the CMB and SNIa data have 
been combined with independent constra in ts from galaxy 
redshif t surveys (e.g. ISpergel et alJ <l2003l) . iTeemark et alJ 
I 2004D ). The cluster X-ray gas mass fraction and galaxy 
redshift survey data break parameter degeneracies in the 
CMB data in similar ways. This is shown in Figure 1101 
which compares the results for a constant equation of state 
model using the Cluster+CMB data (green, dashed con- 
tours) or when combining the CMB data with con s traint s 
from the 2dF galaxy redshift survey jPercival et al.l i2002t) : 
red, solid con t ours. The are the same constraints used by 
ISpergel et all (120031) .) It is clear from the figure that the 
inclusion of the cluster gas fraction data leads to tighter 
constraints on fi m and, especially, wq. However, when the 
SNIa data are also included, the constraints from the clus- 
ters+CMB+SNIa and 2dF+CMB+SNIa data sets are sim- 
ilar, due again to the complementary nature of the vari- 
ous data sets. For the 2dF+CMB+SNIa data, we obtain 
w = -1.02i°:i3 and fi m = 0.304+0.031, which are compa- 
rable to the values in TableLH Note that the 2dF+CMB con- 
straints show n in Fig urc|TU|are somewhat weaker than those 
presented bv lSpergel et al.f (I2003T) . This is due primarily to 
the in clusion of dark energ y perturbations in the present 
study. ISpergel et alJ J2003I) also employ a prior r < 0.3 
which has a small effect on the contours. Such a prior is 
not included here. 

Using our most general dark energy model, we find that 
the combination of 2dF+SNIa+CMB data leads to slightly 
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Figure 11. The 68.3 and 95.4 per cent confidence limits in the (ra s , f2t>/i 2 ) plane from the analyses of SNIa+CMB data (left panel) and 
clustcrs+CMB data (right panel) using the dark energy model with zt = 0.35. Also shown are the thinned samples, where the colours 
correspond to the value of r. Note how the combination of clusters+CMB data alleviates the degeneracies between these parameters. 




Figure 12. The 68.3 and 95.4 per cent confidence limits in the (fi m ,u>o) planes from the analysis of all three data sets combined, 
assuming either a flat prior (solid contours) or including the curvature as an additional free parameter (dashed contours). The left panel 
shows the results for constant w and the right panel for a the Zt = 1 dark energy model. 



weaker and more negative values for wo and larger values for 
w e t, which rule out the cosmological constant at la. How- 
ever, as discussed above, the departure from the cosmologi- 
cal constant becomes insignificant when we also include the 
cluster gas fraction data. 



5 DISCUSSION 

In this paper, we have examined the evolution of dark energy 
using a combination of three cosmological data sets: CMB 
anisotropies, SNIa and the gas fraction in X-ray luminous 
galaxy clusters. Employing a minimum of prior information, 
we have obtained tight constraints on various key parame- 
ters including the dark energy equation of state. Assuming 
that the equation of state remains constant with time, and 
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assuming a flat prior, we measure wo = — 1.05^° 'J2. Employ- 
ing a series of more general dark energy models, we find no 
significant evidence for evolution in the equation of state. 
A cosmological constant is a good description of the cur- 
rent data. For models other than a cosmological constant, 
we have included the effects of perturbations in the dark en- 
ergy component to avoid violating energy-momentum con- 
servation. We have shown that neglecting perturbations can 
lead to spuriously tight constraints on dark energy models, 
especially for the w e t parameter (by up to a factor two for 
our most general model). 

Although each of the data sets used here probes cer- 
tain aspects of dark energy by itself, much tighter con- 
straints are obtained when the data are combined. The 
SNIa and cluster data provide the primary source of in- 
formation on the evolution of dark energy. Of these, the 
SNIa data currently contribute the stronger constraints. 
However, this power can only be utilised once a tight con- 
straint on fi m , in this case provided by the combination 
of the cluster+CMB (or cluster+BBN+HST) data, is in- 
cluded. Using the SNIa data alone, one is hampered by a 
degeneracy betwe en the e quation of state w and_p m e.g. 
iRiess et all J2004D (see also lWang fc Teemarkl ll 2004) ). This 
is the reason that some authors have empl oyed a strong prior 
on f2 m in their studies u sing S NIa data (EI am et aL 
Jassal. Baela. fc PadmanabharJ 120041 : Ijonssori^t^dT 



0.85 




Jain. Alcaniz. fc Devll2004f) ; without such a prior, the SNIa 



data alone cannot provide tight constraints on even constant 
equation of state models, much less models with evolution 
in w. Rather than introducing strong priors, our approach 
has been to use a combination of data sets that are comple- 
mentary in nature and which allow certain key parameter 
degeneracies to be broken. 

One of the ma in parameter degenera cies highlighted 
in previous studies (Corasan iti et al.1 12004'1 is between the 
scalar spectral index n s , the physical baryon density Qbh 2 
and the optical depth to reionisation, r; this degeneracy im- 
pi nges on the measured d ark energy parameters. As noted 



pu rges on the mcasi 
by ICorasaniti et all 



12004ft . the integrated Sachs- Wolfe ef- 
fect in the case of an evolving dark energy equation of state 
increases the importance of this degeneracy with respect to 
constant w models. The left panel of Fig llll shows this de- 
generacy for the case of the SNIa+CMB data, using the 
z t — 0.35 dark energy model. The right panel shows how 
the degeneracy is lessened when the clusters+CMB data are 
used. The combination of clu sters+CMB data a lso leads to 
a tight constraint on Ho (e.g. lAllen et alJ <l2003l) L 

The combination of clusters+CMB data even allows us 
to relax the assumption that the Universe is flat, although we 
note that the computation of MCMC chains in the non-flat 
case is time consuming when one wishes to ensure conver- 
gence with R— 1 < 0.1. (For this reason, we have only carried 
out a limited exploration of non-flat models here.) In order 
to avoid unphysical regions of the parameter space when us- 
ing non-flat models, we have also included a prior on the 
optical depth to reionisation, r < 0.3 , in a similar manner 
toWMAP team JSpergel et al]l200St and ICorasaniti et alJ 
J2004T) . Figure W7\ shows the 68.3 and 95.4 per cent confi- 
dence limits in the (Qm,wo) plane obtained assuming a con- 
stant dark energy equation of state (left panel) and using 
the Zt = 1 model (right panel), with (dashed lines) and 
without (solid lines) the curvature included as a free param- 
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Figure 13. The 68.3 and 95.4 per cent confidence limits in 
the (f2 m ,Q(j e ) plane from the analysis of the combined clus- 
ter+SNIa+CMB data set, with included as a free parame- 
ter. The solid contours show the constraints for constant w. The 
dashed contours show the results for the zt = 1 dark energy 
model. 



eter. Comparison of the curves shows that the uncertain- 
ties in the parameters are increased when the assumption 
of flatness is dropped. However, we still have clear evidence 
that wo < —1/3 and therefore that the Universe is accel- 
erating at late times. For the constant equation of state 
model, we obtain tight constraints on f\ = — 0.0171q'q2i > 
Q m = 0.3Uto°of 6 and Q dc = 0.703t£^. The median' pa- 
rameter values and parameter values at the peaks of the 
marginalised probability distributions (and 68 per cent con- 
fidence intervals) for wo , w ct and f2 m using the non-flat mod- 
els are summarised in Table |2] Figure H"3l shows the results 
in the (f2 m ,S7dc) plane. 

Finally, it is encouraging to recognise the prospects 
for advances in this work over the next 1 — 2 years. Fur- 
ther Chandra observations of X-ray luminous, high-redshift, 
dynamically relaxed clusters should lead to rapid improve- 
ments in the constraints from the X-ray method. Continual 
progress in SNIa studies is expected over the next few years 
and the forthcoming, second release of WMAP data should, 
at the very least, provide an important, overall tightening of 
the constraints. In the long term, the combination of comple- 
mentary constraints from missions such as Constellation-X, 
SNAP and Planck, combining high precision with a tight 
control of systematic uncertainties, offers our best prospect 
for understanding the nature of dark energy. 



ACKNOWLEDGEMENT 

We acknowledge helpful discussions with A. Lewis, S. Bridle 
and P. S. Corasaniti and technical support from R. M. John- 
stone and S. Rankin. The computational analysis was car- 
ried out using the Cambridge X-ray group Linux cluster and 
the UK National Cosmology Supercomputer Center funded 



© 2004 RAS, MNRAS 000.HHTT1 



Constraining dark energy with galaxy clusters, supernovae and the CMB 11 



by SGI, Intel, HEFCE and PPARC. DR is funded by an 
EARA Marie Curie Training Site Fellowship under the con- 
tract HPMT-CT-2000-00132. DR thanks KIPAC, SLAC and 
Stanford University for hospitality during his visit there. 
SWA thanks the Royal Society for support at Cambridge. 
This work was supported by the NASA grant *** and the 
DOE at Stanford and SLAC. JW is supported by the DOE 
and the NASA grant NAG 5-10842 at Fermilab. 



REFERENCES 

Alam U., Sahni V., Starobinsky, A. A., 2004, JCAP, 0406, 
008 

Allen S. W., Schmidt R. W., Ebeling H., Fabian A. C, van 

Speybroeck L., 2004, MNRAS, 353, 457 
Allen S. W., Schmidt R. W., Bridle S. L., 2003, MNRAS, 

346, 593 

Astier P., 2001, Physics Letters B, 500, 8 
Armendariz-Picon C, Mukhanov V., Steinhardt P. J., 

2000, Phys. Rev. Lett. 85, 4438 
Bassett B. A., Corasaniti P. S., Kunz M., 2004, preprint 

(astro-ph/0407364) 
Bean R., Hansen S. FL, Melchiorri A., 2001, Phys. Rev.D, 

64, 103508 

Caldwell R., Dave R., Steinhardt P., 1998, Phys. Rev. Lett., 
80, 1582 

Caldwell R. R., 2002, Physics Letters B, 545, 23 
C aldwell R. R ., Doran M., 2005, preprint 

(astro-ph/0501104) 
Carroll S. M., Hoffman M., Trodden M., 2003, Phys. 

Rev.D, 68, 023509 
Chevallier M., Polarski D., 2001, Int. J. Mod. Phys., D10, 

213 

Corasaniti P. S., Bassett B. A., Ungarelli C, Copeland 

E. J., 2003, Phys. Rev. Lett., 90, 091303 
Corasaniti P. S., Kunz M., Parkinson D., Copeland 

E. J.. Bassett B. A.. 2004, Phys. Rev.D. in press 

(astro-ph/0406608) 
Efstathiou G., 1999, MNRAS, 310, 842 
Ferreira P., Joyce M., 1998, Phys. Rev.D, 58, 023503 
Freedman W., et al., 2001, ApJ, 553, 47 
Gelman A., Rubin D. B., 1992, Statist. Sci., 7, 457 
Hinshaw G., et al., 200 3, ApJ S., 148, 13 5 
Hu W., 2004, preprint (astro-ph/0410680 ) 
Huterer D., Turner M. S., 1999, Phys. Rev.D, 60, 081301 
Jain D., Alcaniz J. S., Abha D., 2004, preprint 

(astro-ph/0409431) 
Jassal H. K., Bagla J. S., Padmanabhan T., 2004, preprint 

(astro-ph/0404378) 
Jonsson J., Goobar A., Amanullah R., Bergstrom L., 2004, 

preprint ( astro-ph/0404468 ) 
Kirkman D., Tytler D., Suzuki N., O'Meara J. M., Lubin 

D., 2003, ApJ S., 149, 1 
Knop R. A., et al., 2003, ApJ, 598, 102 
Kogut A., et al., 2003, ApJ S., 148, 161 
Kuo C. L., et al., 2004, Astrophys. J., 600, 32 
Lewis A., Bridle S., 2002, Phys. Rev.D, 66, 103511 
Lewis A., Challinor A., Lasenby A., 2000, Astrophys. J., 

538, 473 

Linder E. V., 2003, Phys. Rev. Lett., 90, 091301 



Maor I., Brustein R., Steinhardt P. J., 2001, Phys. Rev. 
Lett., 86, 6 

Pearson T. J., et al., 2003, Astrophys. J., 591, 556 

Peebles P., Ratra B., 1988, ApJ, 325, L17 

Percival W.J., et al., 2002, MNRAS, 337, 1068 

Perlmutter S., et al., 1999, ApJ, 517, 565 

Rapetti D., Weller J., 2005, to be published 

Ratra B., Peebles P., 1988, Phys. Rev.D, 37, 3406 

Riess A., et al., 1998, Astron. J, 116, 1009 

Riess A. G., et al., 2004, ApJ, 607, 665 

Saini T. D., Raychaudhury S., Sahni V., Starobinsky A. A., 

2000, Phys. Rev. Lett., 85, 1162 
Spergel D. N., et al., 2003, ApJ S., 148, 175 
Verde L., et al., 2003, ApJ S., 148, 195 
Wang Y., Tegmark M., 2004, Phys. Rev. Lett., 92, 241302 
Tegmark M., et al, 2004, Phys. Rev.D, 69, 103501 
Weller J., Albrecht A., 2001, Phys. Rev. Lett., 86, 1939 
Weller J., Albrecht A., 2002, Phys. Rev.D, 65, 103512 
Weller J., Lewis A. M., 2003, MNRAS, 346, 987 
Wetterich C, 1988, Nucl. Phys., B302, 668 
Zlatev I., Wang L. M., Steinhardt P. J., 1999, Phys. Rev. 

Lett., 82, 896 



© 2004 RAS, MNRAS 000.IT1ITT1 



